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In this paper we implement an improved (error sensitive) Richardson-Lucy deconvolution algo- 
rithm on the measured angular power spectrum from the WMAP 3 year data to determine the 
primordial power spectrum assuming different points in the cosmological parameter space for a flat 
ACDM cosmological model. We also present the preliminary results of the cosmological parameter 
estimation by assuming a free form of the primordial spectrum, for a reasonably large volume of 
the parameter space. The recovered spectrum for a considerably large number of the points in 
the cosmological parameter space has a likelihood far better than a 'best fit' power law spectrum 
up to AXoff ~ —30. We use Discrete Wavelet Transform (DWT) for smoothing the raw recovered 
00 . spectrum from the binned data. The results obtained here reconfirm and sharpen the conclusion 

' drawn from our previous analysis of the WMAP 1st year data. A sharp cut off around the horizon 

I scale and a bump after the horizon scale seem to be a common feature for all of these reconstructed 

OA . primordial spectra. We have shown that although the WMAP 3 year data prefers a lower value of 

matter density for a power law form of the primordial spectrum, for a free form of the spectrum, we 
can get a very good likelihood to the data for higher values of matter density. We have also shown 
that even a flat CDM model, allowing a free form of the primordial spectrum, can give a very high 
likelihood fit to the data. Theoretical interpretation of the results is open to the cosmology commu- 
nity. However, this work provides strong evidence that the data retains discriminatory power in the 
cosmological parameter space even when there is full freedom in choosing the primordial spectrum. 



I. INTRODUCTION 
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. Increasingly accurate measurements of the anisotropy in the temperature of the cosmic microvifave background 
' (CMB) has ushered in an era of precision cosmology. A golden decade of CMB anisotropy measurements by numerous 
experiments was topped by the results from the data obtained by the Wilkinson Microwave Anisotropy Probe (WMAP) 
[l, 2, 3]. Under simple hypotheses for the spectrum of primordial perturbations, exquisite estimates of the cosmological 
' parameters have been obtained from the angular power spectrum measurement by WMAP data combined with 
■ other cosmological observations 0, Q- Precision measurements of anisotropics in the cosmic microwave background, 
[ and also of the clustering of large scale structure, suggest that the primordial density perturbation is dominantly 
OS ■ adiabatic and has a nearly scale invariant spectrum [2, l5| . This is in good agreement with most simple inflationary 
scenarios which predict power law or scale invariant forms of the primordial perturbation 0, 0j @] ■ The data have 
, also been used widely to put constraints on different parametric forms of primordial spectrum, mostly motivated by 
' inflation [ol. [lOl. [ill. [T3.Tl^ . However, despite the strong theoretical appeal and simplicity of a featureless primordial 
t spectrum, it is important to determine the shape of the primordial power spectrum directly from observations with 

minimal theoretical bias. Many model independent searches have been made to look for features in the CMB primordial 
^ power spectrum [Tsl . [TgI [T^ . [isl . [l^ . Accurate measurements of the angular power spectrum over a wide range of 
• multipoles by WMAP has o p ene d up the possibility of deconvolving the primordial power spectrum for a given set of 

] cosmological parameters [lo, IMl l22l . [23l , [23 . [25[ , [29] . Theoretically motivated models that give features in the power 
. spectrum have also been studied and compared in post-WMAP literature [13, IH, [2^, [10, [U, [H, [H, [11] . 
■ - - ' The angular power spectrum, C^, is a convolution of the primordial power spectrum P(A:) generated in the early uni- 
verse with a radiative transport kernel, G(/, fc), that is determined by the current values of the cosmological parameters 
indicated by other cosmological observations. The remarkably precise observations of the angular power spectrum Ci 
by WMAP, and the concordance of cosmological parameters measured from different cosmological observations opens 
up the avenue to directly recover the initial power spectrum of the density perturbation from the observations. The 
error-sensitive Richardson-Lucy (RL) method of deconvolution was shown to be a promising and effective method to 
recover the power spectrum of primordial perturbations from the CMB angular power spectrum (23 |. We have im- 
proved the deconvolution method by factoring out the normalization factors from the iteratively recovered primordial 
spectrum, -P(fc), in the algorithm to remove the artifacts that were present at the two ends of the recovered spectrum 
in our previous work(corrected by template subtraction). For a given set of cosmological parameters, this method 
obtains the primordial power spectrum that 'maximizes' the likelihood to data. 

In this paper we apply the method to the CMB anisotropy spectrum given by WMAP 3 year data. We employ 
Discrete Wavelet Transform (DWT) for smoothing the raw recovered spectrum from the binned data. In this work we 
first present detailed results of an automated computation of the primordial power spectrum for 6 distinct points in the 
cosmological parameter space for flat ACDM models using WMAP 3 year data. Each of these 6 points in the parameter 
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space has specific characteristics of interest. We also present the preliminary results of the cosmological parameter 
estimation optimized over the form of the primordial spectrum in a coarsely sampled volume of the parameter space. 
In this case, instead of simply computing the likelihood for a given model of initial power spectrum, one obtains the 
initial power spectrum that maximizes the likelihood at a point and assigns that likelihood to that point in the space 
of cosmological parameters. However our results for the cosmological parameter estimation, have a coarse resolution 
in spacing of the parameters and is also limited in volume of the parameter space covered. In principle it is possible to 
extend this work to explore the "entire" space of cosmological parameters with high resolution along the lines being 
done routinely. 

In Sections II we review in brief the Richardson-Lucy deconvolution method and the improvements that we have 
made to adapt the method to our problem. In section III we explain the smoothing by discrete wavelet trans- 
form(DWT). The recovered spectrum from WMAP 3 year data for 6 different points in the cosmological parameter 
space is described in Section IV. In Section V we discuss about cosmological parameter estimation by optimizing 
(maximizing) the likelihood over a free form of the primordial spectrum. Finally, in Section VI we discuss our results 
along with concluding remarks. 

II. RICHARDSON-LUCY DECONVOLUTION METHOD 

The Richardson-Lucy (RL) algorithm was developed and is widely used in the context of image reconstruction in 
astronomy [H,!!^. However, the method has also been successfully used in cosmology, to deproject the 3-D correlation 
function and power spectrum from the measured 2-D angular correlation and 2-D power spectrum (stI . [ssj . 

The angular power spectrum, Ci, is a convolution of the initial power spectrum P{k) generated in the early universe 
with a radiative transport kernel, G{l,k), that is determined by the values of the cosmological parameters. In our 
application, we solve the inverse problem of determining the primordial power spectrum, -P(fc), from the measured 
angular power spectrum, Q, using the relation 

G -^G(/,fc,)-P(fc.)- (1) 

i 

In the above equation, the 'target' measured angular power spectrum, Ci = Cf , is the data given by observations, 
and the radiative transport kernel, 

Ak- 

G{l,h)^-^\ATi{h,7jo)\\ (2) 

encodes the response of the present multipoles of the CMB perturbed photon distribution function A^iiki, tjq) to unit 
of power per logarithm interval of wavenumber, k, in the primordial perturbation spectrum. The kernel G{1, k) is 
completely fixed by the cosmological parameters of the 'base' cosmological model. The kernel G{l,k) also includes 
the effect of geometrical projection from the three dimensional wavenumber, fc, to the harmonic multipolc, I on the 
two dimensional sphere. 

Obtaining P{k) from the measured C/, for a given G{1, k), is clearly a deconvolution problem. An important feature 
of the problem is that G^ , G{l,k) and P{k) are all positive definite. We employ an improved RL method to solve 
the inverse problem for P{k) in Eq. The advantage of RL method is that positivity of the recovered P{k) is 
automatically ensured, given G(l,k) is positive definite and C/'s are positive. The RL method, readily derived from 
elementary probability theory of distributions |35|, is an iterative method that can be neatly encoded into a simple 
recurrence relation. The power spectrum P^-^~^^'{k) recovered at the iteration {i + 1) is given by 

p(^+i)(fc)_pW(fc) = pi^)^k)J2Gil,k)^^^^^^ (3) 

I 

where G{l,k) is the normalized kernel (on the I space for all k wavenumbers), Cf is the normalized measured data 
(target) and c|*'' is the angular power spectrum at the i'^ iteration obtained from cl*-* = ^ G(Z, fc) P*^*^ (fc) using 
the recovered power spectrum P^^\k). It is important to remind the reader that due to the formulation in terms 
of conditional probability distributions, the RL method requires the kernel, G(Z, fc), data, G;, and the target vector, 
P(fc), all to be normalized at the beginning. 



3 



^ G = 1; ^^'^w = 1; E = i, (4) 

I k I 

where P^^\k) is the normahzed initial guess model of the primordial spectrum. The normalization factor, = 
G(Z, k) , should be factored out from the final form of the recovered primordial spectrum at the end of the 
procedure, 

P(fc) = p(")(A:)/a, (5) 

where P*^"-'(fc) is the reconstructed primordial power spectrum at the last iteration and P{k) is the final result. The 
Eq. ([3]) and Eq. ([T]) together summarizes the standard RL deconvolution method for obtaining the primordial power 
spectrum P(fc). The final recovered power spectrum is independent of the initial guess p(^^(A:) (see discussion in 
Appendix of ref. (Hi)- 

Due to noise and sample variance, the data Cf is measured within some non-zero error bars u;. The standard 
RL method does not incorporate the error information and tends to iterate to fit features of the noise, as well. In 
our problem, this manifests itself as a non-smooth deconvolved spectrum P(fc) from the binned data that has poor 
likelihood with the full WMAP spectrum data. We devise a novel procedure to make the RL method sensitive to the 
errors tx/ by modifying Eq. ^ to 



p('+i)(fc) -P«(fc) = p(')(/fc) ^ G(/,fc) 



tanh 



Wn2 



(6) 



The idea is to employ a 'convergence' function to progressively weigh down the contribution to the correction p(*+i) — 
p(«) from a multipole bin when cf'' iterate to a value close to CP within the error bar cr/. This innovation significantly 
improves the WMAP likelihood of the deconvolved spectrum. Further improvement of the likelihood to the full Ci 
data is obtained by a subsequent step of wavelet based smoothing described in the next section. In this paper we have 
also improved the IRL method that we used in our previous paper to directly remove the artifacts at the two tails of 
the spectrum. Compared to our earlier implementation, we have eliminated the need to subtract a known model from 
the recovered spectrum to remove the artifacts at the two ends of the spectrum. We have achieved this by a slight 
change in the normalization procedure in the method, that corrects for the effect of very low amplitude of G'(Z, k) at 
very small and very large k. In the previous paper, the normalization factor C^, was hidden in the finally deconvolved 
P(fc) and it was factored out at the end of the process. As we are dealing with an iterative process, every small 
artifact will affect the higher iterations more strongly. To avoid this, we separate the normalization factor C^k from 
the P{k) from the beginning. This small modification gives rise to a big improvement in the method where the final 
form of the recovered spectrum is free from the artifacts at the two ends. Thus we do not need to use any template to 
remove these artifacts (as in our previous paper [22] ) . So in our revised iterative process we have modified Eq. ^ to 



P^'+^\k) - P^'\k) = P«(/fc) V G(;,fc)Cfe-^ tanh^ 

I '-'I 

where the normalization factor, (^k, is explicitly present in the main iterative equation and we do not need to remove 
it from the final form of the P{k) at the end of the process. 



Wn2 



(7) 



P(k) ^ P^^'^k). (8) 

The final form of the recovered spectrum was obtained after smoothing the spectrum by using Discrete Wavelet 
Transform as discussed in the following section. 



III. DISCRETE WAVELET TRANSFORM 



Wavelet transforms provide a powerful tool for the analysis of transient and non-stationary data and is particularly 
useful in picking out characteristic variations at different resolutions or scales. This linear transform separates a 
data set in the form of low-pass or average coefficients, which refiect the average behavior of the data, and wavelet 
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or high-pass coefhcients at different levels, which capture the variations at corresponding scales. As compared to 
Fourier or window Fourier transform, wavelets allow optimal "time-frequency" localization in the real, as well as, 
Fourier domain. The vocabulary of DWT stems from applications in one dimensional time-stream signal trains, but 
has found wide application in signal in other domains and dimensions. Specifically in our case, the 'signal' being 
transformed is the power spectrum, P{k), a one dimensional function of wavenumber, k. 

Wavelets are an orthonormal basis of small waves, with their variations primarily concentrated in a finite region, 
which make them ideal for analyzing localized 'transient' signals. Wavelets can be continuous or discrete. In the 
latter case, the basis elements are strictly finite in size, enabling them to achieve localization, while disentangling 
characteristic variations at different frequencies [39j . This is the preliminary reason for us to employ discrete wavelets 
for our analysis. For more details about DWT and its theoretical basis, see j2^ . 

In this paper, we use DWT to smooth the raw recovered spectrum obtained from the deconvolution using binned 
CMB spectrum data. The raw deconvolved spectrum has spurious oscillations and features arising largely due to the 
k space sampling and binning in I space. The main goal is to reconstruct the primordial spectrum which lead to an 
angular power spectrum with a high likelihood to the entire Ce data at each multipole including the covariance between 
them. The WMAP likelihood of the Ce suffers owing to the spurious oscillations in the spectra on scales smaller than 
£ multipole space bins. The WMAP likelihood improves as the spectra, P{k) is smoothed. We use DWT to smooth 
the recovered spectrum so that WMAP likelihood of the corresponding theoretical Ci is maximized. Strictly speaking, 
the best likelihood obtained in our method is a lower bound leaving open a mathematical possibility of obtaining a 
superior likelihood at the given point of the cosmological parameter space with a different primordial power spectrum 
[5q . Although it is difhcult to establish that the final result is the unique solution with maximum likelihood, but 
numerous variations we have explored does suggest that it is perhaps very close to the best possible result. So we 
can claim that the improved reconstructed likelihoods which we drove for different points in the parameter space by 
assuming a broken scale invariant form of the primordial spectrum, put an upper limits for the best possible results. 

For our smoothing purpose, we use DWT in a systematic way to separate the features of the raw recovered spectrum. 
We map the raw recovered spectrum which has 1400 discrete points to an array of 2^^ points by padding at the two 
ends. By applying the discrete wavelet transform to the input file, we get corresponding 2^^ wavelet coefficients. If 
we apply the inverse wavelet transform to whole set of derived wavelet coefhcients, the orthonormality of the wavelet 
basis will lead back to the raw recovered spectrum. But if we exclude wavelet coefficients above a given level of 
resolution, then the inverse wavelet transform, leads to a smoothed power spectrum compared to raw deconvolved 
spectrum essentially filtering out spurious high resolution structures arising due to numerical effects [57| . 

At the first step, we use the original raw recovered spectrum to calculate the likelihood to the WMAP-3 data. It is 
in fact the inverse wavelet transform of the whole wavelet coefficients which has 2", n = 11 coefficients. We call it the 
recovered spectrum at level 11. In the next step, we cut half of the coefficients and we use only the first 2^° coefficients. 
The derived results (recovered spectrum at level 10) would be smoother than the original spectrum. We continue the 
procedure for all different levels (n=2 to 11), and at each level we calculate the likelihood of the recovered results. The 
recovered result at a level which gives the best likelihood would be our final result. Using discrete wavelet transforms 
has this important advantage of being a well defined smoothing procedure that also retains identify localized features 
in P(k) that contribute a significant improvement to the likelihood. Figure [T] shows the reconstructed P{k) for a 
sample point in the cosmological parameter space, smoothed up to different wavelet levels. The inset of the likelihood 
as a function of smoothing resolution shows that there exists an 'optimal' level of smoothing which 'maximizes' the 
likelihood with respect to data. The blue line in figure [1] is related to the 9th wavelet level where the best likelihood 
to the data is achieved. We reiterate that this allow the entire procedure to be automated so that an "optimal" 
primordial spectrum is recovered given the cosmological parameters. 

IV. PRIMORDIAL POWER SPECTRUM FROM WMAP-3 YEAR DATA 

One of the most challenging questions of modern cosmology is to find an inflationary scenario that satisfies all 
the cosmological observations. The shape of the primordial power of scalar (density) and tensor (gravitational wave) 
generated during inflation are the key observables in this investigation. In our previous paper [22] we introduced 
the Richardson-Lucy deconvolution algorithm to flnd the shape of the primordial power spectrum using the cosmic 
microwave background data for a single point in the cosmological parameter space. In this section we reconstruct 
the primordial power spectrum optimized to get the best likelihood for six different sample points in the cosmological 
parameter space where each of these points has a special importance. We assume a fiat ACDM model and the 
differences between these 6 cases are just in the values of the background cosmological parameters within this sub- 
space of parameters. In the next section we use our automated routine to perform the cosmological parameter 
estimation and explore a coarsely sampled but reasonable sized contiguous volume of the parameter space. In this 
paper, we chose to compare all the reconstructed results with the best fit power law form of the primordial spectrum 
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FIG. 1: Resultant P{k) for a sample point in the cosmological parameter space is shown in the blue curve. The other curves show 
the P(k) recovered at different levels of DWT smoothing. The blue line which is the reconstructed result obtained by retaining 
all wavelet coefficients up to the 9th wavelet level has the best likelihood with A^eff = —18.76 with respect to the best fit 
power-law primordial spectrum in the whole parameter space. We used Ho = 72, 0dm = 0.246, VLh = 0.05, Oa = 0.704, r = 0.06 
as the cosmological parameters. The plot in the inset shows the resultant AXojr of the reconstructed results at different wavelet 
levels. The 'optimality' of the n = 9 level DWT smoothing in this case is clearly demonstrated. 



in the whole parameter space, referred to as "model G" (see Sec. G), rather than comparison with the result from 
the power law form of the spectrum for the same point in the parameter space. This highlights the effect of assuming 
the free form of the primordial spectrum and emphasizes on the significant improvement of the global likelihood. 



A. Cosmological parameters from SDSS 



In this case we consider a flat ACDM cosmological model with cosmological parameters motivated by, and consistent 
with, the results of large scale structure observations from Sloan Digital Sky Survey (SDSS) [41]. We use h — 0.72 
(Hubble parameter), ^dm = 0.246 (dark matter density), fib — 0.050 (baryonic matter density), f^A = 0.704 (A 
energy density) and r — 0.06 (optical depth). These parameters are consistent with the best fit results from SDSS 
with h^ra — 0.213 ± 0.023 and ^hj^m = 0.17 ( where Vim — ^dm + ^b) for a flat ACDM cosmological model. We 
have chosen r = 0.06 throughout the paper, as it is one of the most reliable values for the optical depth at the present 
from observations of Lyman-a forest |42l . l43j . We reconstruct the primordial power spectrum for this point in the 
parameter space. The reconstructed result for P{k) and the resultant Cj"^ and Cj^ are shown in figure [2] (Model 
A). The resultant (including TT and TE polarization spectra) for this point in the parameter space given by the 
reconstructed primordial power spectrum, can improve the effective likelihood by Ax^^j = —18.67 with respect to the 
reference likelihood of model G. 



B. Cosmological parameters from 2DF Galaxy Redshift Survey 

In this case we choose parameters consistent with the results from 2df galaxy redshift survey We use h = 0.63, 
^dm = 0.251, ilfc — 0.041, JIa — 0.708 and t = 0.06. These parameters are consistent with the results from 2df with 
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TABLE I: Different points in the parameter space and the resultant effective likelihood from the reconstructed primordial 
spectrum using WMAP 3 year data. The AXcff is twice the logarithm of the relative likelihood with respect to the best result 
in the whole parameter space by assuming power law form of the primordial spectrum. 



Model 










T 


AXcff 


Model A (compatible with SDSS) 


72.0 


0.246 


0.050 


0.704 


0.06 


-18.76 


Model B (compatible with 2df) 


63.0 


0.251 


0.041 


0.708 


0.06 


-4.38 


Model C (compatible with BAO) 


68.0 


0.229 


0.052 


0.719 


0.06 


-2.93 


Model D (compatible with SN la + BAO) 


72.0 


0.229 


0.046 


0.725 


0.06 


-14.52 


Model E (compare to WMAPl) 


71.0 


0.226 


0.044 


0.730 


0.0 


-13.40 


Model F (compatible with flat CDM) 


50.0 


0.904 


0.096 


0.0 


0.06 


-26.70 



h^lm = 0.168 ±0.016 and n^,/^lrn = 0.185 ±0.046 for a flat ACDM cosmological model. However here we have used a 
marginally bigger value of matter density and marginally lower value of Hubble parameter in compare with the best 
fit result from 2df, but still these parameters are consistent with the 2df constraints within la. The reconstructed 
result for P{k) and the resultant Cf^ and C^^ are shown in figure [2] (Model B). The resultant C/ for this point in 
the parameter space, can improve the effective likelihood by Axis — ~4.38 with respect to the reference likelihood 
of model G. 



C. Cosmological parameters from the results of detection of baryon acoustic peak oscillations 

In this case we consider a fiat ACDM cosmological model with cosmological parameters consistent with the results of 
measurements of the baryon acoustic osciUations(BAO) [i^. We use h = 0.68, ildm = 0.229, Qb — 0.052, Q\ = 0.719 
and T — 0.06. These parameters are consistent with the best fit results from BAO with Vl,nh? = 0.130 ± 0.011 and 
rifeft,^ = 0.024 for a fiat ACDM cosmological model. The reconstructed result for P{k) and the resultant Cj"^ and Cj^ 
are shown in figure [2] (Model C). The resultant Ci for this point in the parameter space given by the reconstructed 
primordial power spectrum, can improve the effective likelihood by A^^g = —2.93 with respect to the reference 
likelihood of model G. 



D. Cosmological parameters from observational constraints on the matter density using SN la and BAO 

data 

Model independent estimation of the matter density by using supernovae [i^, |43] and BAO [1^ data by [i^ is 
used in this section to set the cosmological parameters. The total matter density, independent of the model of dark 
energy is found to be — 0.276 ± 0.023. This value is, in fact, the total sum of dark matter density and baryonic 
matter density. We can also choose the value of baryon density in a way to be consistent with the prediction of big 
bang nucleosynthesis where n^K^ « 0.02 [H, [53|. We use h = 0.72, f^d,„ = 0.229, = 0.046, f^A = 0.725 and 
T = 0.06. These parameters are consistent with the two constraints mentioned above. The resultant C; for this point 
in the parameter space driven by the reconstructed primordial power spectrum, can improve the effective likelihood 
by AXcff = —14.52 with respect to the reference likelihood of model G. The reconstructed result for P{k) and the 
resultant Cj^ and Cj'^ are shown in figure [2] (Model D). 
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E. Cosmological parameters from comparison with the results from WMAP 1 data 

In this case we use the same parameters as we used before in our previous paper [2^ where we introduced the 
improved Richardson-Lucy method and analyzed the WMAP 1 data. This is for comparison between our results from 
WMAP 1 and WMAP 3 years data. As in our previous paper, here also we use h = 0.71, Qdm = 0.226, = 0.044, 
51a — 0.730 and r = 0.0 for the parameters of our background cosmology. The reconstructed primordial power 
spectrum for this point in the parameter space, can improve the effective likelihood by Ax^^j = —13.40 with respect 
to the reference likelihood of model G. This result is consistent with the result by using WMAP 1 data and there is no 
significant difference in features of the reconstructed resultsjl^l. The reconstructed result for P{k) and the resultant 
Cf"^ and Cf^ are shown in figure [2] (Model E). 

F. Cosmological parameters from Standard Cold Dark Matter (SCDM) model 

In this case we assume a Cold Dark Matter universe (CDM) where the energy density of the dark energy is assumed 
to be zero. For h = 0.50, fldm = 0.904, r2{, = 0.096, fi^ = 0.0 and r = 0.06, we could improve the effective likelihood 
by Axcs — —26.70 which shows that by assuming a free form of the primordial spectrum, the standard CDM model 
of a flat universe can be very well fitted to the CMB data alone. Studies by [5l| have shown that a CHDM model 
of the universe (which there is also no dark energy) also can have a good fit to the WMAP 3 years data. It is very 
interesting that for this point in the parameter space we could get a big improvement in the effective likelihood. We 
should note here that for this point in the parameter space, we have set ilbh'^ = 0.024 which is in agreement with big 
bang nucleosynthesis however this point in the parameter space is not well fitted with other cosmological observations, 
like large scale structure observations or the supernovae data. Results are shown in figure [2] (Model F). 

G. Reference model: cosmological parameters from best fit power law to the WMAP 3 year data 

Interestingly, for this point in the parameter space, we could not significantly improve the effective likelihood 
by considering the free form of the primordial spectrum where we used h = 0.732, fldm — 0.1967, fib — 0.0416, 
r^A = 0.7617 and r = 0.089. In fact adding features to the form of primordial spectrum for this point in the 
parameter space could not improve the resultant effective likelihood. We interpret this to arise from the fact the 
cosmological parameters themselves adjust with a large number of degrees of freedom to the fit comparable to the 
freedom in the primordial power spectrum (encoded in a finite number of wavelet coefficients). In short, all the 
cosmological parameters have been chosen in way to give us the best likelihood by strong assumption of power law 
form of the primordial spectrum. In this paper, the resultant XcS °f different models are compared with 
this model as the reference model. 

Different assumed models with their parameters and their resultant likelihoods are shown in table [H 

V. TOWARD COSMOLOGICAL PARAMETER ESTIMATION 

It is very important to note that despite of allowing a free form for the primordial spectrum, not all cosmological 
models (i.e., all points in the parameter space) can be fitted to the data equally well. We clearly show that some points 
in the cosmological parameter space fit the WMAP 3 year CMB data better than the other points, by 'optimizing' the 
likelihood over a free form of the primordial spectrum. We conjecture that the positive definiteness of the primordial 
spectrum does not allow us to fit all the points in the parameter space to the data equally well, and some points will 
have a better fit to the data. In this section we would like to present strong evidence that despite of allowing a free 
form of the spectrum, the derived likelihoods do strongly discriminate between different neighboring points in the 
parameter space. As an example, for a flat ACDM model, we fix the values of Hq = 72, t = 0.06 and f^A = 0.70 and 
we vary the values of fif, and fldm keeping the total fixed at ilm = ilb + i^dm = 0.30 and calculate the likelihood. We 
find a minimum in the value of the XeS around ilf, — 0.050 and fldm = 0.250 which shows that the data prefers this 
combination among the models with fim = 0.30. 

In parameter estimation, other cosmological observations can be used to put constraints on the parameter space. 
For example, some region of the parameter space may be in agreement with CMB data, but being ruled out strongly by 
other observations. However in our example, the best result seems to be well in agreement with all other cosmological 
observations, including large scale structure observations from SDSS [4l|, supernovae data by SNLS 46] and Gold 
[131, detection of baryon oscillations ^4^ and it is also in agreement with big bang nucleosynthesis. In figure [3] we see 
the resultant ^XcS versus different values of baryonic matter (left panel-blue line) and the reconstructed P{k) for the 
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best combinations of fib and i^dm assuming f2m = 0.30 (right panel-blue line). The red curves in figure [3] have the 
same characteristic but f^A — 0.73 and = 0.27. Clearly, 'optimizing' over the primordial power spectrum allows 
us to get significantly higher likehhood (Ax^ = —19.65) for flm = 0.30 compared to = 0.27 (Ax^ = —11.55). 

Here we have only considered two limited 1-D slices in the cosmological parameter space, but our main aim is 
to do the cosmological parameter estimation for the whole volume of the parameter space. However as it has been 
mentioned before, doing a cosmological parameter estimation in the whole parameter space would be computationally 
very expensive. Recently few new methods of parameter estimation have been proposed which are claimed to be 
much faster than the usual methods of Monte Carlo Markov Chain or grid sampling. These new methods may be 
suitable for our purpose but they are still applicable for a one dimensional space and they need to be modified to 
be apphed in our problem [s^. [ssl l&i. [ssj . Though it is still difficult and beyond our abilities to do the cosmological 
parameter estimation in the whole parameter space and with a high resolution, still we can do it for a reasonably 
large contiguous volume of the parameter space. In this part, we present the results for the cosmological parameter 
estimation by fixing the value of r = 0.06 (which we have chosen throughout the paper) and varying the values of 
i^b, ^dm and h in a, large volume of the parameter space. In our analysis we use the following priors: 35 < Hq < 85, 
0.012 < rith^ < 0.030, 0.04 < < 0.30 and < f^A < 1. Any of the given priors for Qb, ^dm and h has been divided 
by eight equally-spaced points wEose combinations will generate our assumed initial sampling in the parameter space. 
CIa is derived from the other parameters assuming a fiat universe. Our motivation to choose these wide priors are 
from the strong limits from other astronomical and cosmological observations. In figure [4] we see the resultant —Ax^g 
(in Z axis and also in color indicated by the tool bar in the upper panel) versus different values of Hubble parameter 
(X axis in both upper and lower panel) and ^Ihh^ (Y axis in both upper and lower panel). The lower panel shows the 
relative values of the JIa m our parameter space (indicated by color in the lower panel) . 

The best likehhood has been derived for Qb = 0.084, fldm = 0.764, 17a = 0.152 and Ho = 50 with Axis = -29.282 
with respect to the reference likelihood of model G. Another point in the parameter space with fib — 0.058, ^Idm = 
0.416, r^A = 0.526 and Hq = 60, which is more compatible with the other cosmological observations, also has a very 
good optimized likelihood with Ax^s — —29.014 with respect to the reference likelihood of model G, which is very 
close to our best result. 

Results also show that for a very wide range of f^A, even when its density is close to zero, we can have a very good 
fit to the WMAP 3 year data which is expected, as CMB data alone are not very much sensitive to the A density. 

Due to computational limitations, we could not do a fine resolution cosmological parameter estimation. Work is 
still in progress and results will be published as soon as we accomplish the computations with the WMAP 5 year 
data. 

We hope that this would provide a completely different perspective to cosmological parameter estimation unbiased 
by prejudices of the early universe. 

VI. CONCLUSION 

We present the reconstruction of an 'optimal' primordial power spectrum for ffat ACDM cosmological models for 
different sample points in the parameter space. 

In section IV we have chosen these sample of points to be consistent with different independent cosmological 
observations, or consistent with special theoretical models. Almost in all cases the recovered spectrums improves the 
resultant effective likelihood significantly in comparison with the best fit power law form of the primordial spectrum 
in the whole parameter space. In section V we generalized our study to a much bigger sample of points where we 
performed the cosmological parameter estimation in a large volume of the parameter space. 

There are some important conclusions that can be drawn from our results. Though the published results from 
WMAP team after release of WMAP 3 year data suggests the lower value of matter density in comparison with 
the other cosmological observations like large scale structure observations from SDSS, supernovae and detection of 
baryon acoustic peak oscillations, our analysis shows that by assuming a free form of the primordial spectrum, CMB 
data can be well fitted to the models with the higher value of matter density compatible with the other cosmological 
observations. Our preliminary studies in the parameter space, show some evidence that by assuming the free form of 
the primordial spectrum models with higher value of matter density are better fitted to the WMAP 3 year data than 
models with low value of matter density. Another important result is that a standard CDM model of the universe, can 
also be very well fitted to the WMAP 3 year data by assuming the free form of the primordial spectrum. In fact for 
a SCDM model, we could get one of our best recovered results by improving the likelihood around Ax'^g = — 27 with 
respect to the best result for the power law form of the spectrum in the whole parameter space (reference likelihood 
of model G). The features of the recovered results for all these points in the parameter space have something in 
common. They all show a sharp cut of around the horizon and some significant features after the horizon scale. It 
has been shown before by {24| that the effect of these features around the horizon are very important in improving 
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the likelihood. These results are in agreement with the results from our previous paper where we studied the WMAP 
1 data. 

Assuming a free form of the primordial spectrum, we derived the best likelihood for a ACDM model at ili, = 0.084, 
^dm 0.764, = 0.152 and Hq — 50 within the explored region with Ax^g = —29.282 with respect to the reference 
likelihood of model G. Another point in the parameter space with fii, = 0.058, fldm — 0.416, Ha — 0.526 and Hq — 60, 
which is more compatible with the other cosmological observations, also has a very good optimized likelihood with 
Axgfj = —29.014 with respect to the reference likelihood of model G. 

The differences between our results in parameter estimation and the results from WMAP team assuming a power 
law form of the primordial spectrum motivates us to work towards precise cosmological parameter estimation allowing 
full freedom to the form of the primordial spectrum, with a higher resolution in the spacing of the parameters and 
also considering the whole parameter space. 

In this work, we use the modestly determined quality of CMB polarization spectra from WMAP (TE and EE) 
simply as a consistency check. However, CMB polarization (EE) spectra from Planck Surveyor is expected to be good 
enough to allow us to extend our deconvolution method simultaneously to CMB temperature and polarization. 
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FIG. 2: Reconstructed primordial spectrum (top panel) and the resultant Cj'^ (middle panel) and Cj^ (lower panel) angular 
power spectra are plotted for 6 different points in the parameter space assuming a flat ACDM cosmological model. Cosmological 
parameters of Model A: Ho = 72, fldm ~ 0.246, Qb ~ 0.05, f^A = 0.704, r = 0.06 and the recovered results for this model gives 
AXoff = -18.76. Cosmological parameters of Model B: Ho = 63, f^d™ = 0.251,^6 = 0.041, = 0.708, r = 0.06 and the 
recovered results for this model gives Axcft = —4.38. Cosmological parameters of Model C: Ho = 68, Qdm ~ 0.229, 57b ~ 
0.052, f^A ~ 0.719, r = 0.06 and the recovered results for this model gives AXcff = —2.93. Cosmological parameters of Model 
D: Ho = 72, Qdm = 0.229,^6 = 0.046, ^a = 0.725, r = 0.06 and the recovered results for this model gives AXeff = -14.52. 
Cosmological parameters of Model E: Ho ~ 71, fld?n ~ 0.226, Qt = 0.044, Qa ~ 0.730, r — 0.0 and the recovered results for this 
model gives AXcB = -13.40. Cosmological parameters of Model F: Ho = 50, Qdm = 0.904, fib = 0.096, Oa = 0.0, r = 0.06 and 
the recovered results for this model gives Ax^s ~ —26.70. Model G is the reference model against which all calculated 
AXcff S are with respect to this model. This represents the best fit power law primordial spectrum in the whole 
parameter space. The red error-bars in the middle and lower panels represents the binned angular power spectrum from 
WMAP 3 year data. The black error-bars at the middle panel at the high £, are from ACBAR experiment [i^. The excess of 
power and the bump in the recovered P{k) at the high k {logk/kh ~), seems to be related to the higher measurements of the 
angular power spectrum at high ^'s in WMAP 3 year data in comparison with the other experiments such as ACBAR. 
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FIG. 3: A 1-D slice {Qm = constant) through the cosmological parameter space demonstrates that the data retains strong 
discriminatory power in the cosmological parameter space even when there is full freedom in choosing the primordial power 
spectrum. Left panel: Plot of AXcff of the reconstructed results with respect to the reference likelihood of model G, by 
assuming free form of the primordial spectrum, for a flat ACDM model with ho — 0.72, r = 0.06 and f^A ~ 0.70 and 
i},n = i}b+Qdm ~ 0.30 for different values of Qb (blue line). The red curve is for similar models except for Q,rn = ^t + i^dm = 0.27. 
Clearly, 'optimizing' over the primordial power spectrum allows us to get significantly higher likelihood (Ax^ ~ —19.65) for 
= 0.30 compared to flm ~ 0.27 (Ax^ ~ —11.55). This demonstrates that even though we allow a free form of the 
primordial spectrum, the data does show very strong preference for particular values of cosmological parameters. Right panel: 
Reconstructed primordial spectrum, P{k), for a flat ACDM model with fij, — 0.050, Sldm = 0.25, ho = 0.72, r — 0.06(blue line). 
For these parameters of Qb and Qdm, we could get the best likelihood for Qm = 0.30. The red line is the reconstructed P{k) for 
a flat ACDM model with ilb = 0.0460, Qdm = 0.224, ho = 0.72, r — 0.06. For these parameters of Qb and Qdm, we could get the 
best likelihood for the Qrn = 0.27. It is clear that the reconstructed P{k) for these two points in the cosmological parameter 
space are very similar. However the resultant AXcff for these two points in the parameter space shows a big difference. 
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FIG. 4; A coarse resolution and limited volume exploration of the cosmological parameter space demonstrates that the data 
retains strong discriminatory power in the cosmological parameter space even when there is full freedom in choosing the 
primordial power spectrum. The resultant — Axfa is shown (in Z axis and also in color indicated by the tool bar in the upper 
panel) versus different values of Hubble parameter (X axis in both upper and lower panel) and Q,bh^ (Y axis in both upper 
and lower panel). The lower panel shows the relative values of the f^A in our parameter space (indicated by color in the lower 
panel). We have assumed here r = 0.06. 



